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Abstract 

A complete scheme for production, cooling, accelera- 
tion, and ring for a 1.5 TeV center of mass muon collider 
is presented, together with parameters for two higher en- 
ergy machines. The schemes starts with the front end of a 
proposed neutrino factory that yields bunch trains of both 
muon signs. Six dimensional cooling in long-period heli- 
cal lattices reduces the longitudinal emittance until it be- 
comes possible to merge the trains into single bunches, one 
of each sign. Further cooling in all dimensions is applied 
to the single bunches in further helical lattices. Final trans- 
verse cooling to the required parameters is achieved in 50 T 
solenoids. 



Table 1 : Parameters of three muon colliders. Note 1 : Depth 
is relative to any nearby low land, e.g. Fox river at FNAL. 
Note 2: Survival is from the end of phase rotation to the 
collider ring. 
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INTRODUCTION 

Muon colliders were first proposed by Budker in 
1969 |1|, and later discussed by others ||3l. A more de- 
tailed study was done for Snowmass 96 |4|, but in none of 
these was a complete scheme defined for the manipulation 
and cooling of the required muons. 

Muon colliders would allow the high energy study of 
point-like collisions of leptons without some of the diffi- 
culties associated with high energy electrons; e.g. the syn- 
chrotron radiation requiring their acceleration to be essen- 
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Figure 1 : (Color) Schematic of the components of a Muon 
Collider. 



tially linear, and as a result, long. Muons can be accel- 
erated in smaller rings and offer other advantages, but they 
are produced only diffusely and they decay rapidly, making 
the detailed design of such machines difficult. In this paper, 
we outline a complete scheme for capture, phase manipula- 
tion and cooling of the muons, every component of which 
has been simulated at some level. 

The work in this paper was performed as part of the 
NFMCC collaboration 0, the recently formed MCTF la, 
and Muons Inc. ITj. 

COLLIDER PARAMETERS 

Table [T] gives parameters for muon colliders at three en- 
ergies. Those at 1.5 TeV correspond to a recent collider 
ring design ||9l . The 4 TeV example is taken from the 96- 
study m . The 8 TeV is an extrapolation assuming higher 
bending fields and more challenging interaction point pa- 
rameters. All three use the same muon intensities and emit- 
tances, although the repetition rates for the higher energy 
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Figure 2: (Color) Transverse vs. longitudinal emittances 
before and after each stage. The nine stages are indicated 
with the numeral 1-9. 

machines are reduced to control neutrino radiation. 

PROPOSED SYSTEM 

Figure [T] shows a schematic of the components of the 
system. Figure[2] shows a plot of the longitudinal and trans- 
verse emittances of the muons as they progress from pro- 
duction to the specified requirements for the colliders. The 
subsystems used to manipulate and cool the beams to meet 
these requirements are indicated by the numerals 1-9 on 
the figures. 

Proton Driver 

The proton driver requirements depend on the muon sur- 
vival estimates that will be discussed in a later section. We 
further assume, from the neutrino factory studies, that pion 
production in the 21 best bunches, at the end of phase ro- 
tation, is 1.7% per proton per GeV. The resulting required 
proton bunches, for different energies, are given in Tb. [2l 
For efficiency in the following phase rotation, an rms bunch 
length of 3 ns is required. The space charge tune shift and 
required longitudinal phase space densities are challenging 
at the lower proton energies, but easier at the higher ener- 
gies. 

Table 2: Proton bunch intensity for three different proton 
energies 

I E(GeV) 8 24 60~ 
Np(10i3) 21 7 2.8 



Production, Phase Rotation, and Initial Cooling 

The muons are generated by the decay of pions produced 
by proton bunches interacting in a mercury jet target O. 
These pions are captured by a 20 T solenoid surrounding 
the target, followed by an adiabatic lowering of the field to 
3 T in a decay channel. 
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Figure 3: Phase spaces during phase rotation a) before 
bunching, b) after bunching, c) after rotation. 



The first manipulation (#1), referred to as phase rota- 
tion fTOl, converts the initial single short muon bunch with 
very large energy spread into a train of bunches with much 
reduced energy spread of which we use only 21. The ini- 
tial bunch is allowed to lengthen and develop a time-energy 
correlation in a 110 m drift. It is then bunched into a train, 
without changing the time-energy correlation, using rf cav- 
ities whose frequency varies with location, falling from 
333 MHz to 234 MHz. Then, by phase and frequency con- 
trol, the rf accelerates the low energy bunches and deceler- 
ates the high energy ones. Figure[3] shows ICOOL [1 IJ sim- 
ulations of the phase spaces before bunching, after bunch- 
ing, and after rotation. 

Muons of both signs are captured and then (#2) cooled 
transversely in a linear channel using LiH absorbers, peri- 
odic alternating 2.8 T solenoids, and pillbox 201 MHz rf 
cavities. All the components up to this point are identi- 
cal to those described in a recent study 1 121 for a Neutrino 
Factory. 
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Figure 4: (Color) One cell of the first 6D cooling lattice. 



6D Cooling Before Merge 

The next stage (#3) cools simultaneously in all 6 di- 
mensions. The lattice |13|, shown in figure (H uses 3 T 
solenoids for focus, weak dipoles (generated by tilting 
the solenoids) to generate bending and dispersion, wedge 
shaped liquid hydrogen filled absorbers where the cooling 
takes place, and 201 MHz rf, to replenish the energy lost in 
the absorbers. The dipole fields cause the lattices to curve, 
forming a slow upward or downward helix (see inset in 
Fig. [5]). The following stage (#4) uses a lattice essentially 
the same as (#3), but with twice the field strength, half the 
geometric dimensions, and 402 instead of 201 MHz rf. Fig- 
ure [5] shows the results of a simulation of both systems us- 
ing ICOOL. Although this simulation was done for circular, 
rather than the helical geometry, it used realistic coil and rf 
geometries. Preliminary studies |[T4ll suggest that the dif- 
ferences introduced by the helical, instead of circular, ge- 
ometries will be negligible. The simulation did not include 
the required matching between the two stages. The simu- 
lations also used fields that, while they satisfied Maxwell's 
equations and had realistic strengths, were not actually cal- 
culated from specified coils. Simulations reported in refer- 
ence |[T3l . using fields from actual coils, gave slightly better 
results. 




turns (c=33 m) 

Figure 5: (Color) ICOOL simulation of 6D cooling in 
stages (#3) & (#4). Inset: long pitch helical geometry of 
(#3). 



Bunch Merge 

Since collider luminosity is proportional to the square of 
the number of muons per bunch, it is important to use rela- 
tively few bunches with many muons per bunch. Capturing 
the initial muon phase space directly into single bunches 
requires low frequency (^30 MHz) rf, and thus low gra- 
dients, resulting in slow initial cooling. It is thus advanta- 
geous to capture initially into multiple bunches at 201 MHz 
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Figure 6: (Color) ID Simulation of merge (#5): a) before 
(blue) and after (red) first rotation, b) after second rotation. 



and merge them after cooling allows them to be recom- 
bined into a single bunch at that frequency. This recombi- 
nation (#5) is done in two stages: a) using a drift followed 
by 201 MHz rf, with harmonics, the individual bunches are 
phase rotated to fill the spaces between bunches and lower 
their energy spread; followed by b) 5 MHz rf, plus harmon- 
ics, interspersed along a long drift to phase rotate the train 
into a single bunch that can be captured using 201 MHz. 
Results of an initial one dimensional simulation of this pro- 
cess is shown in Fig. [S] Work is ongoing on the design 
and simulation of a system with the low frequency rf sepa- 
rated from a following drift in a wiggler system with nega- 
tive momentum compaction to reduce the length and decay 
losses. 



6D Cooling After Merge 

After the bunch merging, the longitudinal emittance of 
the single bunch is now similar to that at the start of cool- 
ing. It can thus be taken through the same, or similar, 
cooling systems as (#3) and (#4): now numbered (#6) and 
(#7). One more (#8) stage of 6D cooling has been designed 
(Fig. [7]), using 10 T magnets, hydrogen wedge absorbers, 
and 805 MHz rf. Its ICOOL simulated performance is 
show in Fig. [H Again, the simulation shown used fields 
that, while they satisfied Maxwell's equations and had real- 
istic strengths, were not actually calculated from specified 
coil configuration. 
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Figure 7: (Color) One cell of the last 6D cooling lattice. 
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Figure 8: (Color) ICOOL simulation of final 6D cooling 
lattice (#8) using 10 T solenoids and 805 MHz rf. 

Final Transverse Cooling in High Field 
Solenoids 

To attain the required final transverse emittance, the 
cooling needs stronger focusing than is achievable in the 
6D cooling lattices used in the earlier stages. It can be ob- 
tained in liquid hydrogen in strong solenoids, if the mo- 
mentum is allowed to fall sufficiently low. But at the 
lower momenta the momentum spread, and thus longitu- 
dinal emittance, rises relatively rapidly. However, as we 
see from Fig.[2l the longitudinal emittance after (#8) is far 




Figure 9: (Color) Results of ICOOL simulations of trans- 
verse cooling in liquid hydrogen in 7 sequential 50 T 
solenoids. 



less than that required, so such a rise is acceptable. Figure|9] 
shows the results of ICOOL simulation of cooling in seven 
50 T solenoids. The simulation did not include the required 
matching and re-accelerations between the solenoids. 

A 45 T hybrid Cu and superconductor solenoid! 15] is 
currently operating at NHFML and an upgrade to 50 T 
is planned, but this magnet uses a lot of power. A 27 T 
solenoid using an 8 T YBCO insert at 4K has recently been 
successfully tested lfT6ll . There is a conceptual design of an 
all superconducting 50 T solenoid ifTTll . 

Acceleration 

Sufficiently rapid acceleration is straightforward in 
linacs and recirculating linear accelerators (RLAs). Lower 
cost solutions might use fixed field alternating gradient 
(FFAG) accelerators, rapidly pulsed magnet synchrotrons, 
and/or hybrid SC and pulsed magnet synchrotrons |T8l. 

Muon Losses 

The estimate of muon losses is very preliminary. The 
simulations assumed Gaussian initial distributions and 
were not very well matched into each lattice, leading to 
larger initial losses. And no tapering of the focus param- 
eters as function of length was included, leading to larger 
losses as the emittances approached their equilibrium. As a 
result, the losses observed were larger than those deduced 
using the performances in the mid range of each simula- 
tion. Table [3] shows the result of an attempt to estimate 
realistic losses, but this remains very preliminary. Since 
it is this estimate that was used to determine the required 
proton driver specifications used above, these too must be 
considered very preliminary. 

Table 3: Calculated transmission tune shifts at different 
stages in the system. 
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Collider Ring 

For the 1.5 TeV c.m.s, a lattice has been developed (ll. 
The parameters as given in Tb. [T] were achieved with ade- 
quate momentum acceptance but with dynamic transverse 
acceptance of only at little over 2cr for the specified final 
emittance. We note however that since luminosity is de- 
pendent on the square of the bunch densities, there would 
be little luminosity loss if the larger amplitudes were colli- 
mated prior to injection into the ring. 



Space Charge Tune Shifts 



CONCLUSION 



For bunches with Gaussian distributions in all dimen- 
sions: 

^ fN^\ P± ave 

where f3i_ ave = {2^^) and = 1.35 IQ-^^ m 
Then at the the ends of a number of stages in this system, 
one obtains the tune shifts given in Tb. IH 



Table 4: Calculated maximum space charge tune shifts at 
different stages in the system. 
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Note that A/"^ is larger at earlier cooling stages to allow 
for losses. The first order shifts can be corrected by increas- 
ing the focus strength, but tune spreads of half the shifts 
cannot be corrected. 

Before the merge, the shifts are small because the num- 
bers of muons per bunch are small. The only 6D cooling 
stage with significant tune shift is the last (#8). Its tune 
accepted Aiy/u ^ 0.7 which is 5 times the calculated max- 
imum full tune spread of ± 7.3%, and is not expected to be 
a problem. 

The tune shifts in the 50 T cooling will be significant 
only during the reaccelerations, where we have assumed 
P±s corresponding to 3 T focusing fields. The design of 
these lattices to accept such tune shifts appears possible, 
although we are clearly nearing the limit. 

ONGOING STUDIES 

There is a serious question as to whether the specified 
gradients of rf cavities operating under vacuum would op- 
erate in the specified magnetic fields. This is under study 
by NFMCC collaboration |5 | and alternative designs using 
high pressure hydrogen gas, or open cell rf with solenoids 
in the irises, are being considered. The bunching and phase 
rotation will be optimized for the muon collider, instead of 
being copied from a neutrino factory. Instead of the slow 
helices, a planar wiggler lattice is being studied that would 
cool both muon signs simultaneously, thus greatly simpli- 
fying the system. The use of more, but lower field (e.g., 
35 T) final cooling solenoids is also under study. Experi- 
ments are underway to demonstrate two of the new tech- 
nologies: mercury target |8|, ionization cooling L19J . Fur- 
ther experimental studies are needed. 



Although much work remains to be done, the scenario 
outlined here appears to be a plausible solution to the prob- 
lems of capturing, manipulating, and cooling muons to the 
specifications for muon colliders with useful luminosities 
and energies, even up to 8 TeV in the center of mass. 
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